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ABSTRACT 



. Radiative torques on irregular dust grains, in addition to producing 

! superthermal rotation, play a direct dynamical role in the alignment of 

interstellar dust with the local magnetic field. The equations governing the 
orientation of spinning, precessing grains are derived; H2 formation torques and 
paramagnetic dissipation are included in the dynamics. Stationary solutions 
(constant alignment angle and spin rate) are found; these solutions may be 
stable ( "attractors" ) or unstable ("repellors"). The equations of motion are 
numerically integrated for three exemplary irregular grain geometries, exposed 
to anisotropic radiation with the spectrum of interstellar starlight. The resulting 
"trajectory maps" are classified as "noncyclic", "semicyclic" , or "cyclic", with 
examples of each given. 

We find that radiative torques result in rapid grain alignment, even in the 
absence of paramagnetic dissipation. It appears that radiative torques due to 
starlight can account for the observed alignment of interstellar grains with the 
Galactic magnetic field. 



Subject headings: ISM: Dust, Extinction - Polarization - Scattering 
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1. Introduction 

The discovery of the polarization of starlight (Hall 1949; Hall & Mikesell 1949; Hiltner 
1949a,b) revealed not only that interstellar dust particles were nonspherical, but that some 
process had brought about large-scale alignment of these grains. The interstellar magnetic 
field was immediately suggested as a grain alignment agent possessing the large-scale 
coherence required by the polarization observations. 

Initial investigations considered magnetic processes for bringing about the grain 
alignment. Both "compass-needle" alignment of ferromagnetic grains (Spitzer & Schatzman 
1949; Spitzer & Tukey 1949, 1951) and paramagnetic dissipation in spinning grains 
(Davis & Greenstein 1949, 1950a,b, 1951) were initially proposed. When HI Zeeman 
splitting (Verschuur 1969) and dispersion and Faraday rotation toward pulsars (Woltjer 
1970) indicated interstellar field strengths of 3-4/xG, it was evident that "compass- 
needle" alignment was negligible, but paramagnetic dissipation on the Davis-Greenstein 
timescale tdg [see eq.(|i~9"|) below] appeared to be viable, particularly if the grains were 
superparamagnetic (Jones & Spitzer 1967). However, study of the statistical mechanics of 
grain alignment (Jones & Spitzer 1967; Purcell & Spitzer 1971) raised questions about the 
ability of the interstellar magnetic field to achieve the observed degree of grain alignment, 
since random gas-grain collisions (and magnetic fluctuations within the grain) would tend 
to oppose the alignment process. 

Martin (1971) pointed out that because of the Rowland effect, charged interstellar 
grains would have a substantial magnetic moment either parallel or antiparallel to the 
angular velocity; this magnetic moment ensured that precession of the grain angular 
momentum around the magnetic field would be sufficiently rapid that the observed 
interstellar grain alignment should always be either parallel or perpendicular to the 
magnetic field direction, regardless of the mechanism responsible for grain alignment. Thus 
the galactic magnetic field could lead to large-scale coherence in observed grain alignment, 
even if the magnetic field itself played no direct role in grain alignment. This conclusion 
was reinforced when Dolginov & Mytrophanov (1976) noted that the Barnett effect implied 
much larger magnetic moments for spinning grains, leading to even more rapid precession 
of J around B . 

A very important development occurred when Purcell (1975, 1979) pointed out that 
grains were subject to systematic torques which, in diffuse clouds, would result in grain 
rotation rates far in excess of those which had been previously assumed to result from 
simple thermal excitation by collisions with gas atoms. Purcell considered three sources 
of systematic torques: H 2 formation at preferred sites on the grain surface, photoelectron 
emission following absorption of a UV photon, and variations in the "accommodation 
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coefficient" over the grain surface; the H 2 formation torques appeared to be most important. 
Purcell observed that if these torques were long-lived, then normal paramagnetic relaxation 
would gradually bring the grain angular momentum into alignment with the magnetic field. 
However, if the grain surface changes on time scales short compared to tdg, substantial 
grain alignment would not occur. Spitzer & McGlynn (1978) analyzed the "crossover" 
process when the grain's rotation rate (in body coordinates) undergoes a reversal due to 
a change in sign of the "Purcell torque", and obtained an estimate for the disalignment 
caused by frequent crossovers. Spitzer & McGlynn concluded that if the active sites for H 2 
formation were short-lived, the disorientation during crossovers was too great for grains to 
be aligned by ordinary paramagnetic dissipation. The crossover phenomenon has recently 
been reanalyzed by Lazarian & Draine (1997), who find that the disalignment per crossover 
can be (paradoxically) substantially suppressed by "Barnett effect fluctuations" , with the 
result that paramagnetic alignment could potentially explain the observed alignment of 
a £ 0.1/zm grains, if no other systematic torques acted to change the grain alignment. 

In a previous paper (Draine & Weingartner 1996a; hereinafter Paper I) we found that 
interstellar starlight exerts very substantial torques on irregular grains. These "radiative 
torques" drive extreme superthermal rotation in grains with the sizes (a ~ 0.2/xm) which 
must be aligned to account for the observed interstellar polarization. In addition, radiative 
torques can act to change the direction of the grain rotation, and must therefore play a role 
in the process of grain alignment. This paper is devoted to a study of this phenomenon. 

The grain geometry and coordinate systems are introduced in §[|. §|3| discusses the 
torques acting on the grain due to the Barnett effect, gas drag, H2 formation, paramagnetic 
dissipation, and absorption and scattering of starlight. 

In §^ we obtain the equations of motion after averaging over precession, and we discuss 
the stationary points and "crossover points" allowed by these equations. "Trajectory maps" 
showing the evolution of the grain orientation and rotational velocity are discussed in §[5], 
where we show that three different classes of grain behavior are possible. Examples of 
each of the three classes of trajectory map are presented. In many cases radiative torques, 
together with magnetic precession of J around B , rapidly bring about alignment of the 
grain angular momentum J with the magnetic field Bo- 

In we summarize the behavior of our three exemplary grains as a function of the 
angle ip between the magnetic field B and the starlight anisotropy direction. Since the 
surface density of active sites for H 2 formation is uncertain, we examine the sensitivity 
of our results to this quantity, and find that, for likely values, radiative torques are more 
important than paramagnetic dissipation for producing grain alignment. In some cases the 
radiative torques do not allow a stationary solution; even in this case the average alignment 
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can be substantial. Our results are summarized in §0. 

2. Grain Dynamics in a Magnetic Field 

2.1. Grain Geometry 

As in Paper I, we consider irregular grains of density p and volume V, with an "effective 
radius" a c g = a_5 cm. In the present paper we report results for three 

specific grain shapes, shown in Fig. [I]. Shape 1, which can be described as an assembly of 
13 cubes, is the target geometry considered in Paper I. Shape 2 can be described as an 
assembly of 11 cubes, with coordinates given in Table 1. Shape 3 can be described as the 
union of 5 overlapping spheres, with coordinates and radii given in Table 2. 

The moment of inertia tensor has eigenvalues I\ > 1% > ^3, with principal axes ai, a2, 
and a3.[] We define dimensionless parameters ay by 

2 

L j = a j^P Va ls ■ (!) 

A sphere has ctj = 1; the irregular grain of Paper I has oi\ = 1.745. In Table || we give ctj 
for shapes 2 and 3. 

The thermal rotation rate for the grain is 

V87rai pa 5 ) V100K/ \ P J 

for rotation around ai with kinetic energy kT/2. 

2.2. Coordinate System 

Consider unidirectional radiation, propagating in direction k. Without loss of generality 
we may choose the unit vectors e 1; e 2 , e 3 defining the "scattering coordinates" such that ei 
is parallel to k, and the magnetic field B lies in the e 1; e 2 plane (see Figures |2|, Let ip 
be the angle between Bo and k; since the sign of Bq does not affect the grain dynamics, we 
need only consider if) e [0, 7r/2]. Let the "alignment angle" ^ G [0, 7r] be the angle between 



1 Axes ai and a.2 are obviously determined only to within a choice of sign. We arbitrarily choose one of 
the two solutions for ai and likewise for a.2, and "freeze" these in body coordinates. The axis a.3 = ai x a.2. 
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the grain principal axis a.i and B . The superthermally- rotating grain will be spinning 
around a^ because of the Barnett moment (see § |3.1| ), the axis a x of the spinning grain will 
rapidly precess in a cone around B , and it is convenient to define a "precession frame" 
using polar coordinates (r, £,0) with Bo as the polar axis. Then (£, 0) corresponds to the 
direction of the grain axis ai: 

ai = e x (cos if) cos £ — sin if) sin £ cos 0) + e 2 (sin if) cos £ + cos if) sin £ cos 0) + e 3 sin £ sin , (3) 

and unit vectors in the directions of increasing £ and are given by 

£ = — ei (sin if) cos £ cos + cos t/> sin £) + e 2 (cos -0 cos £ cos — sin ip sin £) + e 3 cos £ sin (4) 
= ei sin if) sin — e 2 cos if) sin + e 3 cos . (5) 

The orientation angles and <3> in the scattering frame (see Fig. |2|) can be determined from 

0(£, 0) = cos -1 [cos-0 cos£ — sin-0 sin£ cos0] , (6) 



2 tan" 1 



sin O — sin if> cos £ — cos if) sin £ cos ( 



(7) 



sin £ sin 

The grain orientation is also determined by the angle f3 measuring rotation of the grain 
around ai (see Paper I and Fig. 0). In the present study we will assume the grain to spin 
rapidly around a x (see below) and will therefore employ radiative torques obtained after 
averaging over angle /3; henceforth, this angle will not be explicitly mentioned. 



2.3. Grain Dynamics 



When uj 2 3> u) 2 ^, internal dissipation in the grain will tend to bring the principal axis of 
largest moment of inertia a x into alignment with the angular moment J, as this is the state 
of minimum kinetic energy for fixed J. As discussed by Purcell (1979), both viscoelasticity 
and Barnett relaxation are effective. In fact, thermal fluctuations will act to prevent perfect 
alignment of J with a x (Lazarian & Roberge 1997). In the analysis below we will assume 
perfect alignment of ai with J, so that u> = u)ki. Then 



h 



duo - d£ „ d(f> 



r B + r 



rad 



Th 2 



drag 



DG 



where Tb is the torque due to the Barnett moment, T ra d is the radiative torque due to 
starlight, rn 2 is the torque due to H 2 formation, T^rag is the drag torque due to gas atoms 
(and photon emission), and T DG is the "Davis-Greenstein" torque due to paramagnetic 
dissipation. These torques are discussed below. 
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3. Torques 

3.1. Barnett Moment and Precession around B 

Spinning grains develop a magnetic moment fj, antiparallel to u>, primarily due to the 
Barnett effect (Dolginov & Mytrophanov 1976); if the grain is charged there is a small 
correction due to the Rowland effect which we shall neglect. The Barnett magnetic moment 
is 

a* = w , (y) 

g^B 

where V is the grain volume, is the Bohr magneton, g fa 2 is the gyromagnetic ratio, 
and x(0) is the static susceptibility. The magnetic properties of interstellar grains have 
recently been reviewed by Draine (1996); "normal" paramagnetism is expected to result in 
x(0) ~ 10~ 4 , and superparamagnetism would give even larger values of x(0)- 

The torque resulting from the Barnett effect magnetic moment is 

r B = a* x B = — 4>IiQ B u sin£ . (10) 

This torque will cause the grain to precess in the Galactic magnetic field B with a 
precession frequency 

hu 2a ig ii B pa 2 a ~ 5 \ Q,p j V 10 -'/WSJ 

where we have set (7 ~ 2. It is therefore clear that interstellar grains will precess around B 
very rapidly compared to all other timescales except the grain rotation period itself. 



3.2. Gas Drag 

The drag torque is 

T drag = h U T^ ag , (12) 

where the timescale for gas drag is (Paper I) 

na lP a e{i 4 «i / p \ 1/2 (3000 cm^K \ 

where the drag coefficient 5 is of order unity. A sphere has 5 = 1; we estimate 5 fa 2 for 
shapes 1 and 2, and 5 fa 1.5 for the more compact shape 3. 

There is also drag on the grain due to emission of far-infrared radiation; while it 
can be comparable to gas drag for very small grains (Paper I), it is unimportant for the 
a e fj ~ 0.1/im grains we shall be interested in here, so we shall neglect it. 
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3.3. H 2 Formation 

Purcell (1979) estimated the torque due to H 2 formation on the grain surface, after 
averaging over the grain rotation, to be 

Th 2 • an = I (|) / n(H) {2EkTfl 2 al a I p(t) . (14) 

where / 2 is the surface area per H 2 formation site on the grain surface, / is the fraction of the 
arriving H atoms which depart as H 2 , E is the kinetic energy of the departing H 2 molecules, 
and p(t) is a random variable with time averages (p(t)) = 0, (p(t)p(t + r)) = e -r / <L , where 
£l is the "lifetime" of a surface recombination site. We define a characteristic rotation rate 

5 /3X 1 / 3 / / £ \ 1/2 /ra(H)\ I 
^ = 24 UJ U^) (15) 

such that the H 2 torque, averaged over the grain rotation, is 

r H2 = ai ^%(t) . (17) 

T~drag 

Note that wh 2 is proportional to /, the characteristic separation between H 2 formation sites. 
We will usually assume / = 10 A for purposes of illustration, but other values of / will be 
considered below. 



3.4. Paramagnetic Dissipation 

Magnetic dissipation in the grain produces a torque (Davis & Greenstein 1951; Jones 
& Spitzer 1967) 

Tdg = -(tsin^cos^ + ai sin 2 ^)/i^ r D ^ (18) 

where 

^c = ^^ = i.5xiov(^U0^H^V ; (19) 

5K(u)B$ ySgcm^y \ K[u) j \ B ) 

K{lo) = x"(cj)/cj, where x" is the imaginary part of the complex susceptibility. We expect 
K ~ 10~ 13 s for normal paramagnetism and u £ 10 9 s~ 1 (Jones & Spitzer 1967; Draine 
1996). 
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3.5. Radiative Torques 



As in Paper I, we approximate the interstellar radiation field in a diffuse cloud by an 
isotropic component with energy density (1 — j)u rac i plus a unidirectional component with 
energy density 7« ra d- The radiative torque can be written (cf. Paper I) 



2 \ 

Ur add eff A 



(7 [F(£, 0)C + G(£, 0)0 + #(£, 0)ai] + (1 - 7) (^°)ai 



(20) 



F(^,<f>) = (Qr) sin-0 cos£ cos0 — cos?/' sin£) + 

(Qr) •e2(cos , cos£ cos0 — sin£ sin-0) + (Qr) -63 cos£ sin0 (21) 
= (Qr) -ei sin-0 sin — (Q r ) -e 2 cos -0 sin + (Q r ) -e 3 cos (22) 
iJ(£,0) = (Qr) ■e 1 (cos?/> cos£ — sin-0 sin £ cos 0) + 

(Qr)-e 2 (sin-0cos£ + cos -0 sin £ cos0) + (Qr)-e 3 sin£ sin0 (23) 

where the radiative torque efficiency vector Qr(6,$) depends on £, 0, and the angle -0 
between B and k, Qf° is the torque efficiency factor for isotropic radiation, and () denotes 
averaging over the incident radiation spectrum, here assumed to be the interstellar radiation 
field ("ISRF") spectrum (Mezger, Mathis, & Panagia 1982; Mathis, Mezger, & Panagia 

1983) , for which u rad = 8.64 x 10" 13 ergs cm" 3 , A = 1.202/xm. (Qf°)isRF is given in Table 
H for shapes 1-3; we see that it is only ~ 1% of ai • (Qr) isrf(Q = 0), showing that the 
radiative torque due to isotropic starlight is negligible compared to that due to anisotropic 
starlight, unless the anisotropic component is less than ~ 1% of the total background. 
The angles 6(£,0) and $(£,0) are obtained using eq. (|6],|7|). The grains are assumed to be 
composed of material with the dielectric function of "astronomical silicate" (Draine & Lee 

1984) . 

As discussed in Paper I, it is sufficient to compute only Q r (0,O). Expressions for F, 
G, and H in terms of Qr(6,0) are given in the Appendix. Q r (6,0) was computed as 
described in Paper I, assuming $ = and a number of different values of the grain rotation 
angle (3. 

As discussed in Paper I, computation of Qr(O,0) as a function of both and 
wavelength A is very cpu-intensive. Accordingly, we have considered only the three grain 
shapes of Fig. [[], and only a single size (a e ff = 0.2/im) for each shape. 



4. Grain Dynamics 
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4.1. Equations of Motion 

From eq. (B) we obtain three equations of motion: 



dtp 




~dt 


H\lo sin £ 


d£ 




dt 


2I x u 


du 




~dt 


2/1 y 



sin £ cos £ 
^) + (l-7)(Qr°) 



^ sin 2 £ ^HgKQ 



a; 



TDG 



(24) 
(25) 
(26) 



rag 



Because the precession rate [see eq. (|TT|) 1 is so much faster than the other terms in 
eq. (|24|426]) we may immediately average over in eq.( p5 , p6 ): 

di 

dt' 
duo' 

~dF 

where t' = t/r drag , 



7 ^F(£) -(3 sin£ cos £ 



M 



(3 
V 

M 
F(0 



7~drag 

TDG 
^H 2 



1 H(0 + (1 - 7) (Qr°)J + W (*) - (1 + £ sin 2 £ j J 
1 lu/ut , and 



(27) 
(28) 



WradAaeffTaYag 6.83 X 10 4 / U rad 



2IiUt 
2ir Jo 



5 



1 r 27T 

7T 



A 



1/2 



v /im 
1 /" 27r 
2n Jo 



3 gem" 



(29) 
(30) 

(31) 
(32) 



[Note that the parameter (3 in eq. (|29|) is not to be confused with the angle (3 in Fig. |^.] For 
the numerical examples in this paper, we assume conditions appropriate to an interstellar 
diffuse cloud: = 30 cm -3 , n(H)/riH = 1,T — 100K, B = 5/xG. We consider a grain of 
silicate composition with p = 3 gem -3 , K = 10 -13 s, and a_5 = 2. For most examples we 
assume / = 1/3, I = 10 A, and E = 0.2 eV. Thus 

'2' 



(3 
V 
M 



2.87 x 10" 



V2 /2 



272 1 nJ 

ai \ x / 2 



1.49 x 10 5 



1.5 



10A / 

2- 



(33) 
(34) 

(35) 



The functions F(£) and if(£) are plotted in Figures f| and ||for shape 1, a e ff = 0.2/zm, and 
the interstellar radiation field. 
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4.2. Stationary Points 



Because p(t) in eq. ( ]T4] , p!7| , |26| , p8D fluctuates, there are no true steady-state solutions. If 
the surface site lifetime t L 3> r drag , then there will be quasi-steady solutions £ s (p), uj s (p); £ s 
is a zero of the function 



7(p\ - F(P\ /jsin ^ CQS ^ 
Z{i) = m - l+/3sin 2 £ 



7^(o + (i-7)(gn + f 



(36) 



and ijj s is given by 



M 



7#(&) + (l- T )<Qj?°> +VP 



1 + 



(37) 



From Figs. § and [| we see that F and are of comparable magnitude. Since for 



conditions of interest we have 7 « 0.1, M > 1, and /3 < 1 (see eq. |31| , |33D , it is apparent 
from eq. (|36| ) that Z(£) ~ F(£), so that the zeros of Z will nearly coincide with the zeros of 
F . The H2 torques and Davis- Greenstein torques have only a minimal effect on the function 
Z(£), and hence on the angles £ s for stationary solutions (although the angular velocities 
u s are significantly dependent on the H 2 torques, as will be seen below). 

The timescale for approach to alignment angle £ s may be estimated by linearizing eq. 



2g) around (£ s ,u; s 



where 



dt' 
du' 

IF 



A 
B 

c 

D 



A(£-£ a )+B(u/-a/.) 



^MdF 



u' s dt. 

MY u 



(3 cos 2£ s 



7 M^-/3o;>n2£ s 
l + /3sin 2 £ s ) . 



(38) 
(39) 

(40) 
(41) 

(42) 
(43) 



The stationary point (£ 8 ,u} 8 ) is an attractor (i.e., stable) provided both 



A + D < and BC - AD < 



(44) 
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otherwise it is a repellor (i.e., unstable). If stable, there will be two relaxation times for 
approach to the stationary point, given by 

-i -(A + D)±[(A-Dy + ABC} 1/2 

'relax ~ r> • l 40 J 

^'drag 

We take the slower of these two relaxation rates to give the alignment time: 

-(A + D) -U-Dr + iBCF . (46) 



2r d 



rag 



4.3. Crossover Points 

When uj — > 0, eq.(p?D is singular unless F = 0. Therefore uj can only change sign at a 
crossover point £ c where 

HQ = . (47) 

From eq.(^) and eq.(|3^) it is easily seen that £ c = 0, n are always solutions to eq.(f4"7|). The 
crossover point acts as an attractor if 

7 M fdF\ 

M[ 7 ^(e c ) + (i- 7 )(QF°)]+^ \di) > ' ( } 

and as a repellor if eq.(f48D is not satisfied. Physical "crossovers" will only occur at crossover 
points which are attractors. 

We define the "polarity" of a crossover point as positive or negative according to the 
sign of du/dt. From eq. (p8|) we see that 

polarity = sign {m[ 7 #(£ c ) + (1 - 7 )(<2F°>] + Vp} ■ (49) 



5. Trajectory Maps 

5.1. Map Classification 

Provided the rotation is highly superthermal (so that we may assume that ax || J/w), 
and rotation and precession are both rapid compared to the timescale for changes in u or 
£, the grain state is specified by two coordinates: the angle £ between J and B , and the 
rate u of rotation around ai. The dynamical evolution of the grain therefore corresponds 
to a trajectory on the £, u plane. Each point on the £, u plane lies on a single trajectory, 
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except for (1) stationary attractors where many different trajectories have a common end 
point, or (2) crossover attractors, where many different trajectories converge. 

We refer to the set of all trajectories as the "trajectory map" . The £, uj plane naturally 
separates into two regions, uj > and uj < 0, with trajectories connecting the regions only 
at a small number of crossover attractors (a crossover repellor formally admits a single 
trajectory connecting the two regions but this is a set of zero area on the £, uj plane). To 
a considerable extent it is possible to classify trajectory maps according to the types of 
stationary points and crossover points which they possess. Three distinct types of trajectory 
maps are possible. 

5.1.1. Cyclic Maps 

If a trajectory map has no stationary attractors, then the generic trajectory does not 
end. Such trajectory maps are referred to as "cyclic". Fig. [8] is an example of a cyclic map. 

From eq.(|27|) we see that if (5 = (i.e., r DG — > oo: no Davis- Greenstein torque), then 
the sign of d^/dt is independent of uj (except for the sign of uj). Since a closed trajectory 
must have two points with the same £ but opposite signs of d^/dt, it follows that there 
can be no closed trajectories contained within either the uj > or uj < regions. If closed 
trajectories exist, then each such trajectory must pass through u = at an even number of 
crossover attractors, and the trajectory must have equal numbers of crossover attractors of 
positive and negative polarity (usually just one of each). 

For finite tdg closed trajectories which do not cross uj = could, in principle, exist, 
depending on the properties of F(£) and H(£) (although any such trajectories cannot cross 
£ = O).0 However, since (3 -C 1 for interstellar conditions of interest, it seems safe to assume 
that closed trajectories do not occur unless they cross uj = 0. This will henceforth be 
assumed. 



5.1.2. Noncyclic Maps 

If the trajectory map does not have at least one crossover attractor of each parity then, 
following the discussion above, there will be no closed trajectories. We refer to this type of 



2 If a closed trajectory crossed £ = 0, then there would be at least two different points on the trajectory 
with £ = but opposite signs for d£/dt. It is clear from eq.(|27|) that this is not possible since the Davis- 
Greenstein torque term vanishes at £ = 0. 
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map as "noncyclic" . 

In a noncyclic map, every trajectory must terminate at a stationary attractor point, on 
a timescale given by eq. fl46|) ; thus there must be at least one stationary attractor. Figs, fi- 



ll are examples of noncyclic trajectory maps. 



5.1.3. Semicyclic Maps 

The case where there are both closed trajectories but also at least one stationary 
attractor will be referred to as "semicyclic" . A necessary condition for a semicyclic map is 
that there be at least one crossover attractor of each polarity, and at least one stationary 
attractor. If there is no crossover repellor situated between two crossover attractors of 
opposite polarity, then there will be trajectories with \u\ — > connecting the two crossover 
attractors. Therefore a sufficient condition for a semicyclic map is that there be at least 
one stationary attractor, and at least one pair of crossover attractors of opposite parity 
with no crossover repellor in between. Figs. |6| and [?| are examples of semicyclic maps. 

If there is a repellor situated between each pair of opposite-polarity crossover 



attractors, there may be no closed trajectories: Fig. |TT| is an example of such a noncyclic 
trajectory map. We conjecture that whenever a crossover repellor lies between each pair of 
opposite-polarity crossover attractors, the trajectory map will be noncyclic. 



5.2. Examples 

In Figure |6| we show evolutionary trajectories for an a e g = 0.2/im silicate grain with the 
geometry of shape 1 in Fig. [I], for if> = (anisotropic component of the radiation field parallel 
to Bo). We show results assuming p = (no H 2 torques) and tdg = oo (no Davis-Greenstein 
torque): the only torques are due to gas drag and starlight. This is an example of a 
semicyclic map. The arrows in the figure are at intervals of 0.5td ra g — 7.5 x 10 4 yr. We see 
that trajectories beginning at uj > all converge on the crossover point at cos£ = —.873, 
where they enter the region uj < 0. Some of the trajectories in the region uj < evolve to 
the stationary attractor at cos£ = — 1; others lead to the crossover attractor at cos£ = 1, 
where they enter the uo > region. 

As mentioned above, in our analysis the crossover points are singular - we are 
unable to follow a grain's evolution through the crossover event, and therefore we cannot 
predict which of the trajectories emerging from a crossover attractor will be "populated". 
Furthermore, our analysis assumes superthermal rotation with perfect Barnett relaxation 
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so that ;*! || 3/u. This approximation will be accurate when the grain rotation is extremely 
superthermal, {uj /u T ) 2 1, but fails when the rotational kinetic energy approaches thermal 
values. The region {uj/ojt) 2 < 10 is shaded in Figs. |S|-|I(]|; with our present approximations 
we are unable to follow trajectories within this region. 

If we now include H 2 formation torques (with p = +1) and Davis-Greenstein torques 
due to normal paramagnetism (tqg = 1.0 x 10 7 yr) the "flow pattern" changes to that shown 
in Figure |7|. The map remains semicyclic. Figures ^| and [7| are qualitatively similar, showing 
that H2 torques, while by no means negligible, appear to be of secondary importance 
compared to radiative torques for a e fr = 0.2/im grains. 

Figures |6| and [7] each show semicyclic behavior. Some trajectories (e.g., from 
uj = — 300co>r and cos£ = —.9) proceed directly to a stable attractor where they are 
captured. Other trajectories (e.g., from to = SOOut, cos£ = —.8) head for the crossover 
attractor. In Figures ^ and [j] it is clear that even though we are unable to predict how 
the grain will emerge from the crossover point at cos£ = 1, we know that the grain must 
then proceed to the other crossover attractor at cos£ = —.88. At this crossover the details 
matter: some trajectories emerging from the crossover lead to the stable attractor, while 
other trajectories return to the crossover point at cos£ = +1. In principle, it is possible 
that a grain could cycle back and forth between these two crossover points many times; 
it then becomes essential to understand the mechanics (and statistics) of the crossover 
process. Such cyclic behavior will in general be a possibility whenever we have two crossover 
attractors of opposite polarity with no crossover repellor in between. 

Figure § shows the trajectory map for the same grain as in Figure [!], but now for an 
angle ip = 30° between B and the starlight anisotropy. We now have two stationary points 
rather than three, and there are no stationary attractors: the grain must continuously cycle 
between the two crossover attractors, which again have opposite polarity. This is a clear 
example of a "cyclic" trajectory map. 

If we rotate the anisotropy direction to ip — 60° (Fig. |^) the stationary point at 
cos£ = — 1 (which for ip = 30° was a repellor) crosses over into the uj > region and 
becomes a stable attractor. The trajectory map now has a straightforward structure: there 
is only one stable attractor and only one stable crossover. The grain must evolve to the 
uj > region, where it is inevitably trapped by the stable attractor. 

The case ip = 90° (Fig. |T0|) shows new possibilities. Because of the symmetry of the 
problem, this trajectory map is symmetric under reflection through cos£ = 0. There are five 
stationary points: two stationary attractors plus three stationary repellors. There are three 
crossover attractors, at cos£ = ±1 and 0; all have positive polarity. All trajectories lead to 
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a stable attractor with cos 2 £ = 0.514. This is another example of a noncyclic trajectory 
map. 

Finally, to illustrate the diversity of possible behaviors, Fig. |y] shows the trajectory 
map for an a e fr = 0.2/mi silicate grain with the geometry of shape 3 (see Fig. [I]) for & n 
angle ip = between the starlight anisotropy direction and B . For this case there are 
four stationary points, three of which are attractors. There are two crossover attractors 
and three crossover repellors. This is a noncyclic map, but which of the three stationary 
attractors the grain is captured by depends on the initial conditions and (for trajectories 
which pass through the crossover attractor at £ = 130°) on the details of the crossover 
process, which will determine which trajectories emerging from the crossover attractor are 
populated. 



It is by now evident that radiative torques lead to complex grain dynamics, which 
depend delicately on the grain geometry, the angle i/j between the starlight anisotropy and 
the magnetic field B , and on the H 2 formation torque to which the grain is subject. 

In the Rayleigh limit, the linear dichroism of a population of identical grains, each 
spinning around its axis a 1; is proportional to (see, e.g., Lee & Draine 1985) 



where 9 B is the angle between B and the line-of-sight, (R) is the ensemble-averaged 
alignment parameter (or "Rayleigh reduction factor" ) 



6. 



Discussion 



6.1. Grain Alignment 



Ak = (R) sin 



2 8 B n gr 7ral s Q po i(X) 



(50) 




(51) 



and the polarization efficiency factor 



QpoiW = - [Qext(E || a 2 ,k || a 3 ) + Q ext (E || a 3 ,k || a 2 > 
Q ext (~E || ai, k || a 2 ) - Q e xt(E || ai, k || a 3 ) 



(52) 



where Q ex t is the extinction efficiency factor for the grain. While eq. (|50"D for the linear 
dichroism is strictly valid only in the Rayleigh limit (where, of course, Q ex t depends only 
on the direction of E, but not on k) it provides a good approximation even when the grain 
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is not small compared to the wavelength. Hence (R) is a useful statistic to characterize the 
effectiveness of grain alignment. 

Above we have shown 6 examples of trajectory maps, but these are obviously a very 
incomplete sampling of parameter space. To provide a more global picture of the grain 
dynamics, Figures |l^-|T4] show the "terminal" values of the grain rotation rate to and cos£. 
The domains of noncylic, cyclic, and semicyclic solutions are indicated. The region labelled 
"noncyclic?" is where there is at least one stationary attractor, and there are crossover 
attractors of opposite polarity, but where each opposite polarity pair of crossover attractors 
has an intervening crossover repellor. As discussed above, we suspect that such cases are 
noncyclic, but we have not verified this in all cases. 



In Figs. |i4-|T4l the horizontal broken lines show the values of cos£ = ±l/-y/3 = ±0.577 
corresponding to alignment parameter R = 0. Solutions with —0.577 < cos£ < 0.577 are 
"antialigned" , with R < 0, but it is apparent that antialigned stationary solutions are 



relatively rare. For example, Fig. [12] shows that for p = 1 and p = 0, shape 1 is antialigned 
only for ip = 84°; even in this case the antialignment is minimal and there are competing 
stationary solutions with substantial alignment. 



From visual inspection of Figures [12j-4l4], we conclude that in a substantial fraction of 



cases, the grains have noncyclic trajectory maps with stationary solutions with appreciable, 
and often perfect, alignment. 



6.2. Sensitivity to the Magnitude of the H 2 Formation Torque 

The numerical examples considered above assume the characteristic separation between 



H 2 formation sites to be I = 10 A, with the H 2 formation torque given by eq. [14] with 
p = 0, ±1. It is quite possible that the surface density of active sites is lower (/ > 10 A), in 
which case the r.m.s. H 2 torque will be larger. To explore this, we consider shape 3 with 
a = 0.2/xm. In Figure 15 we show the stationary solutions for pi = ±30,50, 100 A. For 
\pl\ £ 20 A the H 2 formation torque dominates the grain spinup, and we no longer have 
cyclic trajectories (all crossovers have the same polarity). However, we still have multiple 
stationary points in many cases. The values of R for these cases do not appear to be 
very sensitive to the value of pi. Furthermore, even though the H 2 formation torques may 
dominate radiative torques in driving superthermal rotation, radiative torques may still be 
more important than Davis- Greenstein torques in bringing about alignment of the grain 



angular momentum with Bq. This is evident in Figure 16, showing the alignment time 



as a function of I for ip = 60°. Since the Davis-Greenstein torque is proportional to u, in 
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the limit of large I the alignment time approaches the Davis-Greenstein alignment time. 
However, if the active sites for H 2 formation are characterized by I £j 300 A (i.e., > 1 active 
site per 10 5 A ) then radiative torques apparently dominate the grain alignment process. 



6.3. Cyclic Maps 

When the trajectory map shows cyclic behavior, the degree of grain alignment will 
depend on what fraction of the time the grain spends with different values of £. Suppose we 
have two crossover attractors, at £ c i and £ C 2, with polarities + and — , respectively. Consider 
a trajectory beginning near £ cl , at £ cl + <5£ and u = 8u. The total time to reach the next 
crossover is 

At = [ c2 diidi/dty 1 . (53) 
JU+se, 

The fraction dP of the time spent in an interval d cos £ is given by 

dP 1 (54) 



c/cos£ At sin^\d^/dt\ 

In Figure [H] we show the distribution computed for 6 trajectories for grain shape 1 and 
ip = 30° (see Fig. ^). The upper panel shows four trajectories emerging from the crossover 
attractor at cos£ = 1; the lower panel is for two trajectories emerging from the crossover 
attractor at cos£ = — 1. For each trajectory we compute the time-averaged alignment factor 

r 1 dP 

(R)= / dcos£- -R(0 . (55) 

J-i acos£ 

We see that quite different values of (R) are found for the different trajectories, but in 
this case at least there is a tendency for the average alignment (R) to be positive. Thus, 
although we are not yet able to predict which trajectory (or distribution of trajectories) will 
be taken by grains emerging from the crossover attractors, there is at least an indication 
that even cyclic behavior of the kind shown in Fig. || may result in significant grain 
alignment. 



6.4. Future Work 

The present paper is an exploratory study of the role of starlight torques in the 
dynamics and alignment of interstellar grains. Three grain geometries have been studied, 
but for only one size, a e ff = 0.2/im. The important role of starlight torques has been 
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demonstrated, and we conclude that a ~ 0.2yum grains can often be effectively aligned with 
the magnetic field B by this process. 

A number of issues remain to be investigated: 

f . It is important to examine the dependence of this mechanism on the grain size. In 
Paper I we showed that radiative torques are relatively unimportant for a e s ^ 0.05yum 
grains with the geometry of shape 1; based on those results it seems likely that 
radiative torques will not be able to produce alignment of a O.lfim grains, but 
future work should examine the dependence on grain size in greater detail, with 
the important goal of understanding the observed minimal alignment of grains with 
a < 0.1/im (Kim & Martin 1995). 

2. The present study has been restricted to the dielectric function of silicate material. 
It will be of interest to determine whether other grain materials, e.g., graphite, will 
behave similarly. 

3. We have not attempted to examine the grain dynamics when the rate of rotation 
approaches "thermal" values during the "crossover" process. A detailed study of the 
crossover process is essential for understanding the fate of grains having cyclic or 
semicyclic trajectory maps. 

These issues will be addressed in future investigations. 

7. Summary 

Our principal results are as follows: 

1. As anticipated in Paper I, anisotropic starlight incident on an interstellar grain 
produces a torque which, in addition to spinning the grain up to superthermal rotation 
rates, directly affects the orientation of the grain relative to both the anisotropy 
direction and the direction of the local magnetic field B . Accordingly, previous 
studies of grain alignment, which neglected these torques, are incomplete. 

2. We obtain equations of motion for the grain which are valid in the limit where the 
grain rotation period and the precession period are both short compared to the 
timescale for the grain to change its angular velocity and orientation relative to the 
magnetic field. 
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3. The equations of motion have been integrated for three different grains and many 
different initial conditions. For a given grain and a given angle if) between the radiation 
anisotropy direction and B , we produce a "trajectory map" which shows how the 
grain will evolve under the combined effects of starlight, H2 formation torques, gas 
drag, and paramagnetic dissipation (see Figs. |B|-p]). 

4. Under some conditions the trajectory maps contain a "stationary attractor"; 
trajectories coming near this point are permanenty captured. The stationary attractor 
often (but not always) corresponds to a state of perfect alignment of the grain angular 
momentum J with the magnetic field B . 

5. The trajectory map may contain one or more "crossover attractors" where (formally) 
the grain rotation may reverse. Unfortunately, the assumptions underlying our 
dynamical study break down near such crossover points, and the present study is 
unable to describe the dynamical evolution of the grain near the crossover point. 

6. Trajectory maps can be classified as "noncyclic", "cyclic", and "semicyclic" , based 
upon whether or not there are stationary attractors, and upon the locations and 
polarities of crossover attractors, and the locations of crossover repellors. 

7. In diffuse clouds, for grains with effective radii a eff £ 0.1/im, these starlight torques 
dominate the grain dynamics, if the separation between active sites for H 2 formation 
is / ^ 20 A. Even when H 2 torques dominate the grain spinup (/ £ 20 A) the radiation 
torques dominate the grain alignment process for I £ 300 A. 

8. Based on study of three grain geometries, we conclude that radiative torques due to 
anisotropic starlight appear to be the primary mechanism responsible for the observed 
alignment of interstellar dust grains. 

We are grateful to R.H. Lupton for the availability of the SM plotting package, and to A. 
Lazarian and L. Spitzer, Jr., for helpful discussions. This research was supported in part by 
NSF grant AST-9219283 to BTD, and by an NSF Fellowship to JCW. 

A. Appendix: The Functions F(£,<j>), </>), and H {£,<))) 

The expressions (pl|-p3D for F(£,<f>), G(£,<j>), and H(£,<fi) are given in terms of the 
torque efficiency vectors (Qr(0,$)). As discussed in Paper I (eq. 41), it is possible to 
obtain Qr(©, 3>) from Qr(©, 0), so that the time-consuming calculations required to obtain 
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Qr can be restricted to $ = 0. Thus, we obtain G(£, 0) and $(£, 0) from eq.(^|,0), and then 
obtain F(£,</>), G(£,<j>) and H (£,</>) from Q r (6,0): 



+ 



(Qr(e,o 
(Qr(e,o 
+ (Qr(e,o 
= (Qr(e,o 
+ (Q r (e,o 
+ (Qr(e,o 
= (Qr(e,o 
+ (Qr(e,o 



•ei [— sin0 cos£ cos0 — cos^ sin£] 

• e.2 [cos $ (cos ip cos £ cos — sin ip sin £ ) + sin $ cos £ sin 0] 

•£3 [cos $ cos £ sin0 + sin $(sin sin £ — cos ip cos £ cos 0)] (Al) 

• ei [sin^ sin0] 

• e 2 [sin $ cos — cos $ cos ip sin 0] 

• 63 [cos $ cos + sin $ cos ^ sin 0] ( A2) 

• ei [cos-0 cos£ — sin-0 sin£ cos0] 

■e.2 [cos$(sin?/>cos£ + cos ip sin £ cos0) + sin <£> sin £ sin 0] (A3) 
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Table 1. Coordinates of Constituent Blocks for Shape 2 
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Table 2. Coordinates of Constituent Spheres for Shape 3 
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Table 3. 


Principal Axes a^ 


and factors ocj 


for Shapes 2 and 3 
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Table 4. Radiative Torque Efficiency Factors for Shapes 1-3, a eff = 0.2/im 



shape (Q^isrf ai • (Q t )isrf(® = 0) 



1 

2 
3 



1.99 x 10" 4 
-1.10 x 10~ 4 
2.38 x 10~ 4 



7.01 x 10~ 2 
-2.46 x 10- 2 
-1.39 x 10~ 2 
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Fig. 1. — Three representative grain geometries. Shape 1 is the geometry studied in Paper 
I. 
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Fig. 2. — Grain orientation showing "scattering coordinates", O, <3>, and (3 characterizing 
the grain orientation relative to the incident radiation, k is the direction of propagation of 
the incident radiation. 
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Fig. 3. — Grain orientation showing "alignment coordinates" £ and 0. B is the static 
magnetic field, and ip is the angle between B and k. 
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Fig. 4. — The precession-averaged function affecting grain alignment due to radiative 
torques, where £ is the angle between the grain rotation axis ai and the magnetic field B . 




Fig. 5. — The precession-averaged function H(£) for grain with shape 1 and a c g = 0.2/xm. 
H(£) describes the spinup torque associated with anisotropy of the radiation field, where £ 
is the angle between the grain rotation axis ai and the magnetic field B . 
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Fig. 6. — Trajectories on the cos£, u plane, for ip — 0, with only radiative torques and 
gas drag acting. The grain has the geometry of shape 1 in Fig. [I] and an effective radius 
a e fj = 0.2yum. The gas is assumed to have nn = 30 cm -3 and T = 100K. The radiation 
field is assumed to have u ra d/nukT = 2.09, and an anisotropy factor 7 = 0.1. Arrows are at 
intervals of 0.5td rag = 7.5 x 10 4 yr. There are three stationary points: two repellors (crosses) 
and one attractor (circle). There are two "attractor" crossover points at cos£ = —.873 and 
1 (the crossover point at cos£ = — 1 is a repellor). As discussed in the text, the analysis is 
uncertain in the region (uj/ujt) 2 < 10, which has been shaded. 
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Fig. 7. — Same as Fig.|| but now including torques due to H 2 formation (with / = 1/3, 
I — 10 A, p — 1, and n(H)/nn = 1), and paramagnetic dissipation with K = 10~ 13 s, 
and B = 5/xG. There are three stationary points: one attractor (circle) and two repellors 
(crosses). 
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Fig. 8. — Same as Fig. [7], but for if) = 30°. For this case there are no attractors, and two 
repellors. The grain trajectory presumably ends up cycling around and around, crossing the 
u) = axis at the "crossover" points at £ = and 180°, both of which are "attractors". 



Labels A-F designate trajectories examined in Fig.|17 
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Fig. 9. — Same as Fig. |7|, but for if) = 60°. For this case there is one attractor and one 
repellor. 
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Fig. 10. — Same as Fig. 0, but for ip = 90°. For this case there are two attractors (at 
cos£ = ±0.71 and three repellors. 
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Fig. 11. — Trajectory map for shape 3 and ip = 0, an example of a map which has crossover 
attractors of both polarities, but with no closed trajectories: all trajectories terminate on 
one of the three stationary attractors. The three crossover repellors are at the points labelled 
byR. 
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Fig. 12. — ujIujt and cos£ at attractor stationary points for an a cS = 0.2yum silicate grain 
with shape 1 in a diffuse cloud, as a function of the angle ip between the radiation anisotropy 
and the magnetic field B . A starlight anisotropy 7 = 0.1 is assumed. Results are shown 
for H 2 torques characterized by / = 1/3, I — 10 A, and p = +1,0, —1 [see eq. 0l4j)]. In the 
shaded region there are no attractors, and all trajectories are "cyclic" . Open symbols indicate 
that the trajectory map is "semicyclic" (see text), solutions with cos£ > 0.577 or < —0.577 
are "aligned", with R > 0; solutions with —0.577 < cos£ < 0.577 are "antialigned" , with 
R < 0. 




Fig. 13. — Same as Fig. IT2] but for shape 2. 
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shape 3, a pff =0.2/xm 
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Fig. 15.- 



Same as fig. 



14] but for H 2 formation torques with p = 1 and I = 30, 50, 100 A. 
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Fig. 16. — The alignment time for a c g = 0.2/im grain with shape 3 of Fig. [L], for angle 
ip = 60° between B and the starlight anisotropy, as a function of the surface length scale I, 
where the surface density of active H2 formation sites is l~ 2 . The alignment time is shown 
for H 2 formation torque given by eq. [14], directed parallel (p = +1) or antiparallel (p = — 1) 



to the grain axis ai. We see that for I < 10 3 A, the alignment time is significantly smaller 
than the timescale tdg for alignment by paramagnetic dissipation. For I < 50 A, grain 
alignment takes place in £ 10 6 yr. Note that for I ^ 400 A two alignment times are shown, 
corresponding to £ = and 180°; for £ = the radiative torques assist alignment, while for 
£ = 180° paramagnetic alignment must overcome the radiative torques. 
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Fig. 17.— Plot of the fractional time distribution over cos£ for selected trajectories 
beginning near one crossover point and ending at the other, for an a c g = 0.2/im grain 
with shape 1 and ip = 30° (see Fig. |8|). Distributions are shown for 6 trajectories indicated 
by letters A-F in Fig.[8]. For each trajectory the value of (R) is given. 



